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ABSTRACT 

A  method  has  been  devised  for  integrating  the  pressure-scaled  water- 
vapor  and  C0„  absorber  masses  from  the  surface  to  the  a  -  levels  used 
for  prediction  in  the  Naval  Postgraduate  School  primitive  equation  model. 
Using  empirical  expressions  for  atmospheric  absorptivity,  scattering- 
reflectivity,  cloud-reflectivity  and  earth-surface  reflectivity,  the 
useful  solar  insolations  absorbed  at  earth  and  in  the  key  atmospheric 
layers  above  the  earth  have  been  formulated.   The  terrestrial  cooling 
effect  at  earth  and  in  these  same  key  atmospheric  layers  have  also  been 
formulated  using  recently  published  emissivities  for  the  joint  effects 
of  water  vapor  and  C0_.   As  a  result,  the  radiative  heating  (cooling)  at 
key  levels  may  be  ascertained  for  input  into  the  atmospheric  thermodyna- 
mic equation.   In  addition,  application  of  atmospheric  boundary-layer 
modelling  permits  determination  of  the  surface-layer  turbulent  transports 
(in  the  vertical)  of  sensible  and  latent  heat  at  the  earth's  surface. 
These  turbulent  transports  act  in  such  a  way  as  to  provide  a  heat  balance 
at  the  earth's  surface. 
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1.  Introduction 

The  radiation  package  described  here  is  primarily  an  outgrowth  of  that  used 
in  the  UCLA  general  circulation  model  (Gates  et  al.  1971),  although  some  inno- 
vations from  the  corresponding  model  of  NCAR  (Oliger  et  al.  1970)  and  of  the 
Geophysical  Fluid  Dynamical  Laboratory  (Manabe  &  Strickler  1964)  have  also 
been  introduced. 

The  description  here  of  the  radiation  package  and  its  bearing  upon  the  heat- 
budget  of  the  earth  and  atmosphere  is  similar  to  that  in  use  by  Fleet  Numerical 
Weather  Central  (Kesel  &  Winninghoff  1972)  in  application  to  the  FNWC  primitive- 
equation  (P.E.)  forecast  model.   The  most  obvious  similarities  are  in  the  identi- 
cal use  of  the  same  five  Q-levels  (with  0  =  p/TT)  for  purposes  of  data-input,  and 
also  of  similar  criteria  for  the  occurrence  of  a  one-layer  cloud  deck.   The  main 
difference  in  the  present  from  that  of  the  FNWC  model  lies  in  the  inclusion  of 
water-vapor  and  CO  reduced  optical  masses  in  the  Q-layer  0.4  to  0.2.  This  in- 
clusion  is  made  possible  by  a  realistic  method  of  extrapolation  of  the  water- 
vapor  specific  humidity  to  levels  above  Q=0.4(see  Appendix  A).   A  more  clearly 
defined  system  of  integrations  for  the  computation  of  net  flux  of  terrestrial 
radiation  is  employed  in  this  study  (compared  to  that  of  FNWC)  by  means  of  new 
empirical  absorptivities  of  water  vapor  and  CO   (Sasamori,  1968). 

The  vertical-structure  parameters  needed  to  specify  the  properties  of  the 
radiating-absorbing  model  atmosphere  sounding  are  given  in  Section  2.   The 
solar-radiation  package  is  described  in  terms  of  the  model  in  Section  3,  and 
the  terrestrial  radiation  package  is  developed  in  Section  4. 

2.  The  radiative  vertical-structure  parameters, 
(a)   Computation  of  absorber  masses. 

The  five-layer  model  of  Fig.  1  is  similar  to  that  in  use  in  the  FNWC 
primitive  equation  model.   Since  the  half-layer  a  -  levels,  denoted  by  the 
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Fig.  1.   Definition  of  the  water  vapor  and  CO-  absorber  mass  parameters  in 
the  model.   Moisture  and  temperature  parameters  are  introduced  at  the  odd- 
levels  10.0a  as  in  the  analysis  scheme.   Levels  are  identified  by  their 
values  on  the  10. OQ  scale,  while  layers  are  identified  by  their  level 
boundary  indices  [e.g.  (8,10)],  in  parentheses. 


dashed  lines  are  the  primary  levels  for  data-analysis  of  W,  T,  and  q 
(=  mixing  ratio),  there  are  a  total  of  10  levels  important  for  radia- 
tional  computations  in  the  model.   This  happens  because  for  each  of 
the  five  levels  where  q  is  known,  it  is  possible  to  determine  the 
pressure-scaled  water  vapor  absorbed  mass  AU(2k-l)  and  (similarly  that 
of  CO  ,  AC(2k-l ) ) , associated  with  the  center  of  each  of  the  indicated 
five  odd  numbered  levels.   The  appropriate  formulas  are 

q2k-l  Ap   P2k-1   72 
AU(2k-l)  =  ~^-± (_£iL_L)-/z  (2_i) 

8        Po 
for  water-vapor  absorber  mass  in  the  layer  2k-2  to  2k.   Here  Ap  =  .2TT 
is  a  constant  of  the  analysis  scheme,  and  q„,  ,  is  the  mean  value  of  the 

water  vapor  mixing  ratio  centered  in  the  layer  (2k-2,  2k).   The  pressure- 

0.72 
ratio  factor  (p9,   /p  )     is  the  proposed  scaling  factor  for  collisional 

line-broadening,  after  Moller  (1964).   This  formulation  due  to  Moller, 

disposes  of  the  necessity  of  temperature  scaling  in  the  usual  Lorentz 

line-broadening  formula  (c.f.,  Danard,  1969),  where  the  pressure-ratio 

scaling  is  linear  and  a  temperature-ratio  scaling  is  accomplished  by 

T  -.5 
(— )  *   times  the  mass  within  a  layer. 

o 

For  C0_,  Moller  also  proposes  pressure-scaling  similar  in  form  to  (2-1) 

but  with  a  pressure-ratio  exponent  of  0.65.   For  simplicity  a  single 

pressure  exponent  of  .72  will  be  adopted  for  both  water-vapor  and  C0« 


ab 


sorber  masses.   Thus  the  layer-reduced  absorber  mass  of  CO-  in  N.T.P 


2   i 
cm  per  cm  column  is 

-4   An       p2k-l   °*72 
AC(2k-l)    =  3.14  xlO4^     (^iL-i)  (2-2) 

gPo  po 


In  both  Eqs.  2-1,  2-2  the  subscript  "o"  denotes  an  N.T.P  condition,  e.g., 

p   =  1013.25  mb,  and  p   is  the  standard  density  at  pressure  p  =  p   and 

T  =  273.16°K. 
o 

For  computation  of  radiative  transfer,  the  total  absorber  mass  (pro- 
perly pressure-scaled)  above  a  designated  pressure  level  is  needed.   For 
our  purpose  here,  it  is  most  convenient  to  sum  all  absorber  masses  rela- 
tive to  the  earth's  surface,  where  2k  =  10.   Thus  the  reduced  absorber 
masses  of  Eqs.  2-1,  2-2  when  applied  to  the  layer  (8,  10)  become 

q9     P9   72 

^)^s,io-3-i4gr°"4Ap^>72 

o         o 
The  dual  notation  of  (2-3)  suggests  the  further  procedure  to  be  followed 
for  incrementing  absorber  mass  over  specified  layers  in  the  atmosphere. 
The  notation  of  (2-4) indicates  the  procedure  for  incrementation  of  mass 
adopted 

U6,10  "  U8,10  +  U6,8 

C6,10  =  C8,10  +  C6,8  <2"4> 

Here  the  method  involves  incrementing  absorber  mass  only  upward  over 
successive  layers  in  the  five-layer  atmosphere.   The  full  atmospheric  depth 
of  absorber  mass  is  likewise  denoted  by  the  symbolism 


U0,10  ==  f  _  U2k-2,2k 
k=5 


o 


C0,10  ~  ^=5  C2k-2,2k  (2_5) 


where  each  term  U   2  ok  is  similar  in  form  to  U_    ,  with  the  subscript 
9  replaced  by  2k-l. 

Since  it  is  understood  that  the  integration  of  absorber-mass  proceeds 
only  from  the  surface  upward,  the  second  subscript  "10"  in  (2-5)  will  be 
deletedc   In  line  with  these  concepts,  the  following  table  summarizes  the 
parameters  (and  their  elevations)  needed  at  each  gridpoint  to  perform  a 
radiation  calculation. 


Table  1.   Notation  for  the  distribution  of  radiative 
properties  at  even  levels  over  gridpoints  in  the  five- 
layer  model. 
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Since  the  temperature  and  humidity  analyses  are  carried  along  only  at 
the  odd  levels  a  =.9,  .7,  .5,  .3,  .1,  it  is  generally  necessary  to  extra- 
polate a  temperature  at  the  surface, T- ~,  using 


10 


Tg  =  .5(3T9  -  T7) 


(2-6) 


The  temperatures  at  the  other  even-subscripted  levels  are  available  by 
linear  interpolation  between  temperatures  at  odd-levels  of  10a.   In  general, 
no  effective  temperature  at  level  "0"  is  available,  and  in  fact  the  topmost 
available  temperature  analysis  gives  T, ,  which  will  in  general  be  taken 


This  will  continue  to  be  the  case  unless  stratospheric  temperature-  and 
moisture-modelling  are  undertaken  in  future.   Stratospheric  modelling 
should  be  highly  correlated  with  tropospheric  features  below. 


as  the  effective  temperature  of  the  topmost  absorber-mass  increments  of 

water  vapor  and  CO-  (Uq  and  C„),  respectively. 

2(b)   Vertical  extrapolation  of  the  water-vapor  mass  to  upper  levels. 

Smith  (1966)  has  shown  the  existence  of  a  strong  correlation  between 
q(p)  and  the  surface  value  q1f.  through  the  extrapolative  formula 

fllEl  =  (£_  )*■  (2-7 

q10     p10 

Here  X  is  determined  from  the  water-vapor  sounding.   In  connection  with 

the  present  study,  the  Northern  Hemisphere  soundings  tabulated  by  London 

2 
(1957)  have  been  tested   to  obtain  a  best-fit  technique  between  In  q  and 

In  p  so  that  an  equation  of  form 

a&l  =  (£-)x  (2-8) 

q5      P5 

might  be  obtained.   The  values  of  q(p)  corresponding  to  a  =  0.1,  .3,  .5, 
.7,  .9  were  used  to  solve  for  the  profile-parameter  X,  which  varied  with 
latitude  and  season.   The  correlation  coefficient  which  related  In  q  to 
In  p  was  determined  with  such  high  precision  that  use  of  (2-8)  for  extra- 
polation to  a  =  0.1  was  considered  feasible  when  q(p)  was  available  only 
at  the  levels  p  =  0.3,  .5,  .7,  .9. 

With  q(l)  known,  the  scaled  absorber-mass  U_  modeled  after  U_    of 

2.  o ,  10 

Eq.  2-3  was  formulated  for  use  in  the  model.   In  some  cases  values  of  q„ 
will  also  be  missing  but  such  values  are  subject  to  retrieval  from  Eq.  2-8 
using  the  mixing  ratio  values  from  Q^  0.5. 

The  procedure  of  water-vapor  extrapolation  to  at  least  0^.3  should 
be  considered  a  part  of  the  model  presented  here  with  the  relevant  analysis 


2 

By  LT.  R.  J.  Plante  in  his  1972-73  Master's  thesis. 


routine  attached  as  a  subprogram.   The  values  of  U~  are  quite  small,  but 

the  transmissivities  in  the  terrestrial  spectrum  are  appreciable  in  any 

case. 

2(c)   Criterion  for  the  existence  of  clouds  at  gridpoints. 

In  both  the  solar  and  terrestrial  radiative  transfer  problems,  it  is 

essential  to  know  the  fractional  cloud-cover  CL  at  each  gridpoint.   This 

model,  like  that  of  FNWC,  permits  only  a  single  cloud  layer  to  form.   The 

criterion  for  the  amount  CL  of  cloud-cover  in  existence  at  any  gridpoint 

is  given  by  Smagorinsky ' s  (1960)  algorithm  (applicable  at  approximately 

700  mb): 

CL  =  2.0  (  J  )a=  ?   -  0.7  (2-9) 

s 

Thus  at  every  gridpoint  there  is  a  cloud-covered  area  CL  of  fractional 
amount  0  <  CL  ^  1.0,  and  a  clear-area  fraction  0^  1  -  CL  ^  1.0. 

When  the  fractional  cloud-cover  CL  exists,  the  cloud  is  assumed  to 
extend  everywhere  within  the  grid  area  square  between  O  =  ,8  to  o  =  0.4. 

The  occurrence  of  CL  >  0  of  cloud  amount  profoundly  affects  the  treat- 
ment of  both  the  solar  (Sec.  3   )  and  terrestrial  (Sec.  4   )  radiative 
transfer  problems. 

An  improved  cloud  model  permitting  the  introduction  of  a  two-layer 
cloud  system  both  of  which  fulfill  Smagorinsky-type  requirements  is  en- 
visaged for  future  model  development.   Such  improvements  would  follow  a 
procedure  somewhat  similar  to  that  of  Oliger  et  al.  (1970).   This  refined 
radiative  cloud  treatment  would  permit  CL..  of  intermediate  cloud  between 
0=.8  and  a=.6,  and  CL  of  high  cloud  between  0=.4  and  a=.2„   In  the  model 
of  Oliger  et  al.,  the  high-cloud  layer  is  ascribed  a  solar  transmissivity 
(or  alternatively,  a  reflectivity)  intermediate  between  moist  air  and  a 


mid-level  cloud  deck.   Moreover,  the  distribution  of  the  upper  level 

cloud  elements  is  assumed  to  overlap  the  lower-level  ones  in  a  random 

fashion. 

2(d)   Integrated  absorber-mass  formulas 

If  the  result  of  (2-2)  is  summed  from  p(2k)  =  (,2k)p  to  p  =  p 

using  straightforward  integration,  one  obtains 

p=p 

C(2k)  =  3.14  x  10"4^  V     (E-)-72^ 

8po  J      Po     Po 
°(.2k)p 


C(2k)  -  3'1^21°"  H[l  -  (Pf^)1'72]   (W 3)  (2-10) 

po 


where  H  =  7.995  x  10  cm  is  the  height  of  the  homogeneous  atmosphere. 
Since  the  surface  pressure  TT  may  differ  considerably  from  p  ,  it  is 
convenient  to  adapt  Eq.  2-10  to  the  a  -  levels  of  Fig.  1.   This  is 
easily  done  by  setting  the  ratio 


£^  =  U  0(2k) 

p      p 
*o     *o 

Thus   C(2k)   becomes 

C(2k)  -  3aV7f4H  (J  )U72   [1  -  a,.1"72]  (2-U) 

po 

where   a  =  0.0,    0.2,..., 0.8,    1.0  spans   the  5-layers   shown   in  Fig.    1. 

For  water  vapor,    a   similar   integration  of  Eq.    (2-1)   over  the   layer 

(10,2k)    gives   the  reduced  absorber  mass    in  this   layer: 

u(2k)  .  (2  )1-72(T-^i)  V\n4  [a^1-72  -  a        !•"]  (2-12) 

o  °     n=5 

In  (2-12)  q2n  i  ^s  tne  variable  mean  mixing  ratio  of  water  vapor  in 

1  72 
the  layer  (2n-2,2n).   The  parameter  (TT/p  )    '  is  constant  at  a  given 

gridpoint. 


The  computation  of  C(2k)  and  U(2k)  then  depends  upon  the  values  of 

1.72 
CL,    '  at  level  2k  at  each  gridpoint.   The  following  tabular  values 

1.72 
of  a  *    are  useful  in  these  computations 


Level 

2k 

0 

2 

4 

6 

8 

10 

Q2k 

0.0 

.2 

.4 

.6 

.8 

1.0 

°2k 

0.0 

.06277 

.20680 

.41535 

.68126 

1.0 

With  U(2k)  and  C(2k)  given  by  (2-12,  2-11),  these  mass  parameters  increase 
from  the  earth's  surface  upward. 


3.   The  solar-radiative  heating  package. 

(a)  Partition  of  the  solar  spectrum  at  earth. 

Following  Joseph  (1966),  a  solar  constant  of  2.00  ly  min   is  assumed, 
but  is  subject  to  4$  attenuation  by  oxygen  and  ozone  absorption  above  the 
tropopause.   Thus,  an  effective  solar  constant 

S  =  1.92  ly  min" 
has  been  assumed  entering  the  level  O   =  0.2. 

The  effective  insolation  on  any  given  Julian  date  D  at  level  2,  Fig.  1, 

is  then  obtained  from 

r   "2 
F  =  S(£  )   cos  9  (3-1) 

m 

Here  r/r  is  listed  as  a  function  of  the  date  D  in  Smithsonian  Table  169 
m 

(List,  1963).   In  addition,  the  zenith  angle  6  is  available  for  any  given 

time  step  (GMT),  at  latitude  Cp  and  longitude  X,  using  a  subroutine  of  the 

FNWC  radiation  package  (Kesel  and  Winninghoff,  1972,  Appendix  B). 

Joseph  (1966)  has  devised  a  computationally  simple  partition  of  the 

solar  insolation  into  the  part 

FA  =  0.349  F. 
A 

subject  to  water-vapor  absorption  but  not  to  Rayleigh  scattering  in  the 
atmosphere.   The  remaining  insolation  of  the  partition,  denoted  as 

F_  =  0.651  F 
is  subject  only  to  Rayleigh  scattering,  but  not  to  absorption  in  a  clear, 
moist  atmosphere.   The  presence  of  cloud-decks  introduce  cloud-reflecti- 
vities, however,  in  both  the  F  and  F   solar-energy  partitions.   In  any 
case,  the  F  -energy  comprises  all  wavelengths  X  ^  0.9  (J  m  since  here  water- 
vapor  absorption  is  negligible,  while  the  F  -category  comprises  solar 


10 


wavelengths  X  ^  0.9  (J.  m  where  absorption  by  water  vapor  and  CO  bands 
in  the  near  infrared  is  the  dominant  attenuation  process  acting  on  the 
solar  beam.   In  summary,  the  partition  just  devised  gives  rise  to  the 
parts 


r  "2 

FA  =  0.349S  (-  )   cos 
A  r 

m 

r   -2 
F_  =  0.651S  (-  )   cos 

m 


(3-2) 
(3-3) 


The  component  (3-2)  contributes  to  atmospheric  heating  in  any  vertical 
column  whenever  cos  Q  >  0.   The  component  (3-3)  contributes  only  to  earth- 
surface  insolation,  after  scattering  and  reflective  attenuations  in  the 
column  above  earth  are  considered  (and  when  cos  6  >  0). 

In  some  climatological  and  general  circulation  investigations,  e.g. 
Manabe  and  Strickler  (1964),  it  has  been  found  useful  to  compute  the 
time-mean  cos  0  for  use  in  (3-2,  3-3)  for  the  particular  date  under  con- 
sideration.  The  following  table,  after  Manabe  and  Moller  (1961)  lists 

3        n 
daily  mean  values  of  cos  6  as  a  function  of  latitude  and  time  of  the  year. 

In  addition,  Table  2  lists  the  duration  of  the  day  in  hours  for  these 

same  arguments. 


TABLE  2.   The  seasonal  and  latitudinal  distributions  of  the 
length  of  the  daytime  given  in  parts  of  24  hours  and  those  of 
the  weighted  mean  values  of  cos  6. 


°Lat. 

Fractional  1 

ength  daytime 

Weighted  mean 

cos  0 

Apr. 

July 

Oct. 

Jan. 

Apr. 

July 

Oct. 

Jan. 

5 

.508 

.517 

.500 

.496 

.625 

.587 

.614 

.591 

15 

.521 

.537 

.492 

.471 

.618 

.601 

.579 

.549 

25 

.533 

.562 

.483 

.450 

.599 

.593 

.524 

.474 

35 

.546 

.596 

.471 

.421 

.558 

.567 

.458 

.393 

45 

.562 

.637 

.454 

.362 

.501 

.521 

.379 

.317 

55 

.596 

.708 

.437 

.321 

.423 

.453 

.282 

.203 

65 

.629 

.837 

.404 

.208 

.345 

.369 

.176 

.106 

75 

.750 

1.000 

.329 



.241 

.311 

.071 



85 

1.000 

1.000 





.168 

.318 





sec  0  is  very  nearly  equal  to  the  reciprocal  of  cos  0. 
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3(b)   Disposition  of  F  -insolation 

A.   Clear-sky  F  -insolation 

Joseph  (1966),  found  that  the  Rayleigh-scattering  reflection  values 

at  sea  level  (after  Coulson,  1959)  could  be  fitted  by  least  squares. 

The  result  was  a  best-fit  regression  of  form 

a  (tt,6)  =  .085  +  .25074  log  (-  sec  6)  (3-4) 

Po 

where  Tt   is  the  surface  pressure  and  p   =  1013.25  mb.   The  value  of  (X  is 

o  s 

constrained  never  to  exceed  an  upper  limit  of  0.99  even  for  large  sec  6. 

3(b)B.   Earth-surface  albedo 

Another  reflective-type  parameter  to  be  considered  is  that  of  surface 

albedo  (X  ,  assumed  constant  over  bare  land  as  oc  =0.14  (except  over 
g  g 

deserts') ,  but  variable  with  respect  to  cos  0  over  the  ocean: 

a  (SEA)  =  MAX  {.06,  .06  +  .54(.7  -  cos  0)|  (3-5) 

The  value  of  Qt  in  the  Sahara-Arabian  desert  is  also  considered  a  function 
g 


of  cos  0  through 


a  (DES)  =  .8  -  .6  cos  0  (3-6) 

g 


The  gridpoints  which  lie  in  these  deserts  are  easily  located  in  the  FNWC 

model  as  corresponding  to  I  ^  45  and  latitude  Cp  <  arc  sin  .57.   Also 

within  the  northern  snowline  (determined  from  monthly  climatological 

surface  temperatures  T..  ^  -  3  C),  a  is  obtained  from  the  empirical 

ij  g 

result  (Posey  and  Clapp,  1964): 

(SNOW)  =  0.4  H  1.0  +  — ==== 


I      (CLAT  "  5-y 
ag(SN0W)  =  0.4  |l.O  +  — ^ -V  (3-7) 


(CLAT-  45'>  +  <CLAT  ■  5'>' 
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Here  C    is  the  number  of  latitude  degrees  poleward  of  the  mean  northern 

snowline.   A  similar  expression  applies  for  the  south-polar  ice-snowline. 

These  values  of  CL  have  been  drawn  from  either  the  Kesel-Winninghoff  (1972) 

or  Gates  et  al  (1971)  documentations,  respectively. 

At  any  clear-sky  gridpoint,  the  F_- insolation  is  subject  to  the  Ray- 

leigh  sky  reflectivity  0.     given  by  (3.4),  and  the  appropriate  earth-surface 

s 

reflectivity  OL  .   Considering  the  possibility  of  multiple  reflections 

O 

between  earth  and  atmosphere,  each  of  which  diverts  downward  the  fraction 

a,  of  the  earth-surface  reflectance  Cl  ,  the  retained  insolation  in  the 
s  g' 

earth  is  I   given  by 

so        J 

I  (10)  =  F  (1-a  )[i  +  a  a  +  . . . (a  a  )n+  . . .] (l  -  a  ) 

so         s     s        s  g         s  g  g 

or 

F  (1  -  a  )(l  -  a  ) 

Iso(10>  =  S  1  -aa L  <3"8> 

s  g 
3(b)C.   Cloudy-sky  insolation  in  the  F_-region 

(i)   In  the  cloud  layer 

We  make  use  of  the  Rayleigh  reflection  coefficient  a  of  (3-4)  applied 

only  over  air  above  the  cloud-top  at  p  =  .4TT.   Since  the  cloud  layer  is 

composed  of  large  water  droplets,  it  is  subject  primarily  to  Mie-type  or 

diffuse  back-scatter  which  gives  rise  to  the  cloud-reflectivity  parameter 

R_  (approximately,  R_  =  0.5),  and  also  to  the  downward  scattering  which 

ultimately  must  emerge  through  cloud-base  at  p  =  0o8TT. 

The  two  sky  reflectivity  factors  a  (- sec  6)  due  to  air-scattering 

o 
and  Rp  due  to  cloud-reflectance,  go  on  independently  with  some  downward 

reflection  following  each  cloud-to-sky  reflectance.  Hence,  an  effective 

air-cloud  transmissivity  through  the  cloud-top  similar  to  that  implied  by 
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(3-8)  applies  for  F  -insolation  penetrating  the  cloud-top: 

¥Jl  ~  ac    ^  sec   e>Hl  -  RJ 
s  s        p  L 

F_(4)i  =  2 (3-9a) 

S   L 

Since  the  cloud-penetrating  insolation  F  r(4)i  is  not  subject  to  absorp- 
tion  within  the  cloud,  but  only  to  forward  Mie-scattering,  the  fraction 
emerging  as  transmitted  through  the  base  is  likely  to  have  the  mean  slant- 
path  for  diffuse  radiation,  so  that  F   becomes 

Sv 


V1  -  a«<Hr  sec  6n)][1  -  RJ 

S       S   p        D         L 

F„r(8)i 2 (3-9b) 

Here  sec  6  =  1.66  has  replaced  sec  0  as  in  most  cases  of  diffuse-sky 
slant  paths  (Katayama,  1966).   Wherever  ol  appears  either  in  the  numerator 
or  in  the  denominator,  the  value  of  TT  in  (3-4)  has  been  replaced  by  .4TT 
and  sec  6  by  sec  6  =  1.66.  The  value  F   of  (3-9b)  then  gives  the  trans- 
mitted  diffuse  beam  entering  the  subcloud  layer.   In  this  layer,  F  I 

Sv 

will  be  subject  to  further  multiple  reflections  to  be  discussed  below. 
3(b)C.  (ii)   Partition  of  F  r(8)l  in  the  subcloud  layer. 

S  \j 

The  assumption  is  made  here  that  Rayleigh  scattering  by  air  molecules 

occurs  in  the  air-layer  between  p  =  0.8TT  to  p  =  1.0TT.  Furthermore,  the 

mass  of  air  involved  in  this  layer  is  0.2TT,  and  the  angle  6  is  still 

subject  to  the  diffuse  slant-path  condition  arc  sec  0  =  5/3.  As  a  result, 

the  subcloud  scatter-ref lectivity  a  becomes,  by  analogy  with  Eq.  3-4: 

s 

a  (8-10)  =  .085  +  .25074  log  (r^    |)  (3-10) 

s  p    j 
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The  diffuse  beam  F  c(8)i  which  first  enters  the  subcloud  layer  is 

subjected  to  a  sequence  of  reflections  depicted  schematically  in  Fig.  2. 

Consider  the  beam  B^,  associated  with  the  first  entry  of  F  _(8)i  at 

point  B..  of  the  subcloud  layer  (Fig.  2).   The  fraction  a  is  turned  up- 
1  s 

ward  during  the  first  transit  and  the  fraction  (1  -  a  )  reaches  the 

s 

earth  at  G,  .   In  addition  to  a  ,  the  further  fraction  a  (1  -  a  )  is  also 
1  s  gv     s 

turned  upward  from  the  first  entering  beam  B.G..  before  earth-absorption 

occurs.   The  beam  B,G.  therefore  has  been  decomposed  into  the  parts 

(1  -  a  )(1  -  a  )  F  „       into  earth  at  G, 
s       %'      sC  1 

[a  +  a  (1  -  a  )]F  _       starts  along  G1 -B„ 
s    g      s   sC  12 

The  latter  part  is  assumed  to  reach  the  point  B_,  where  the  subsequent 
multiplicative  fraction  R  starts  down  from  B  towards  G_.   The  inter- 
vening subcloud  scattering  layer  reduces  the  transmission  along  B~G9  by 
the  factor  (1  -  a  );  hence  at  G~,  we  have  the  downcoming  second  diffuse 
beam  insolation  contribution 

[a  +  a  (1  -  a  )]  •  (l  -  a  )  r  F_(8)i 
s    g      s  s   L  sC 

with  a  further  decomposition  into  the  parts 

[a  +a(l-a)](l-a)Ra  f  _  —  reflected  at  G. 
L  s    gv     s        s        C  g  sC  2 

[a  +  a  (1  -  a  )](1  -  a  )(1-  a  )  R^F  _ absorbed  at  G„ 

sgs        sgCsC  2. 

Continuation  of  this  process  of  successive  reflections  between  ground 

and  cloudbase  yields  the  final  cumulative  absorption  at  the  earth,  denoted 


as  I  r(10) 
sC 


IsC(10)  =  d-ccs)(l-ag)FsC+[as4<ig(i-as)](l-as)(i-ag)RcFsC 
+  [as-KXg(i-as)](l-as)3(l-ag)R2  ag  FsC+ +...., 
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and  on  through  higher  terms.  Finally,  upon  simplification,  one  obtains 

FsCi^s+  V1-0^  (1-as)(1-VRC 

In  the  above  accounting,  certain  small  atmospheric  reflection  terms 
from  the  subcloud  layer  to  earth  or  to  cloud  have  been  ignored  as  small 
compared  to  the  air-transmission  factor  (1-a  ). 

It  is  to  be  stressed  that  I  p(10)  given  by  Eq.  3-11  is  the  absorbed 

insolation  at  earth,  resulting  from  the  F  -insolation  of  Eqs.  3-3  and  3-9b. 

The  parameter  a  always  arises  as  a  modification  of  Eq.  3-4  depending 

upon  the  depth  of  the  scattering  layer.  (X  and  R  are  subject  to  vari- 

8      ^ 

ability  depending  upon  the  surface  and  cloud  albedos,  although  R  -  0.5 

is  assumed  as  a  first  approximation.   In  the  FQ-wavelengths  no  atmospheric 

absorption  occurs.   The  weighted  mean  earth-surface  absorption  then 

becomes 

I  (10)  =  (1  -  CL)I   +  CL  I   (10)  (3-12) 

s  so       sC 

where  I   is  given  in  (3-8),  and  I   (10)  in  (3-11). 

SO  S  \j 

A  related  subject  of  some  interest  concerns  the  leakage  back  to  space 
of  the  subcloud  insolation  (Eq.  3-9b)  after  upward  cloud- transmittance. 
Cloud- leakage.   It  is  of  interest  to  compute  the  solar-leakage  through 
the  cloud  base  subject  to  the  approximation  noted  in  the  paragraph  just 
below  Eq.  3-11.   Reference  to  Fig.  2  indicates  that  the  sum  total  beam- 
transmittances  at  the  cloud-base  involve  solar  energy  leaks  of  sequential 
amounts  indicated  by  the  following  sum  of  terms: 

FsC,BT   =  WV  V1-*."   1(1"RC}   +  RCag(1-as)2(1   "  V 

+(Rca  )2(l-as)4(l-Rc)+...| 
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or 

F-  RT=  Fq  [a  +  a  (1-a  )](1-R )[ - r]  (3-13) 

sC»B    sC  s   g    s     C     a  2 

C  g    s 

Comparison  of  (3-13)  with  (3-11)  indicates  that  the  leakage  is  closer 
to  the  order  of  magnitude  of  the  second  (smaller)  term  of  (3-11).   Choos- 
ing a  w  0.05,  and  average  values  a  =  0.14,  R  =  0.5,  indicates  the  leakage 
s  8  ^ 

to  be  about  0.1  F  _(8)l  or  about  3$  of  the  extraterrestrial  insolation  of 

8  Li 

Eq.  3-1.  This  leakage  represents  an  additional  component  to  be  added  to 
the  planetary  albedo  in  cloud-covered  areas  only.  Again  it  is  to  be 
stressed  that  none  of  the  F  radiation  is  absorbed  by  the  atmosphere  in 
this  model. 

3.(c)   Disposition  of  F  -insolation 
(A)  Clear  areas 

Here  absorption  of  the  F  -component  of  the  insolation  is  considered. 

4 
Use  of  the  Manabe-Mdller  solar  absorptivity  function  a(u  sec  6)  in  the 

form 

a  =  .0946  (U2  2k  sec  9),JUJ  (3-14) 

will  be  used  extensively  in  this  discussion  as  well  as  in  that  of  the 

cloudy-sky  case  discussed  below.  The  form  of  Eq.  3-14  indicates  that  the 

absorption  function  a  is  applied  over  the  corrected  water  vapor  path 

between  levels  2  and  2k  (see  Fig.  1)  along  the  zenith-path  of  slant  angle  8. 

Actually  (3-14)  is  applicable  to  the  entire  solar  energy  F.     Restricting 

absorption  to  the  ratio  0.349F  requires  that  the  multiplicative  constant 

4 
The  symbol  a,  denotes  the  fractional  absorptivity  function,  while  A(2,2k) 

denotes  the  absorbed  energy  between  levels  2  and  2k. 
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of  (3-14)  be  adjusted  to  (.0946/0.349)  =  .271.  Hence  the  absorbed  insola- 
tion in  the  layer  (2,6)  becomes  simply 

A(2,6)  =  .271  (U2  6  sec  0)-3O3FA  (3-15) 

Similarly,  the  absorbed  insolation  in  the  layer  (2,10)  is 

A(2,10)=  .271(U2  1Q  sec  0)*3O3FA  (3-16) 

The  absorbed  insolation  in  the  layer  (6,10)  is  then  given  by 

A(6,10)=  A(2,10)  -  A(2,6)  (3-17) 

Finally,  the  difference  formed  by  subtracting  A(2,10)  from  F.  gives  the 

A 

direct  transmission  of  F  -insolation  to  the  earth's  surface.   This  direct- 

A 

transmission  [F  -  A(2,10)]  is  denoted  by  F  (10)  and  is  partially  reflected 

(the  fraction  CI  is  applicable).   The  remainder  (1-CX  )F  (10)  is  absorbed  at 
g  g   A 

earth.   The  F.-insolation  absorbed  at  earth  is  given  by 

IA£10)=  FA[l-.271(U2>10sec  9) * JUJ] (l-ttg)  (3-18) 

In  the  foregoing  derivation,  solar  absorption  was  computed  only  below 
level  2  of  Fig.  1,  as  a  4$  stratospheric  reduction  of  the  solar  beam  was 
previously  applied  in  the  layer  above  level  2.  Likewise  no  accounting  of  the 
terrestrial  spectral  absorption  in  this  topmost  layer  will  be  made.  This 
implies  the  assumption  that  the  layer  a^  .2  is  in  continual  radiative  heat 
balance  at  all  gridpoints. 
3(c). B.   F  -insolation;  Cloudy-sky  areas 

In  the  layer  (2,4),  the  Manabe-Moller  absorption  is  employed  in  a 
manner  similar  to  Eq.  3-15  apart  from  the  differing  water-vapor  mass  UL  , : 
A(2,4)=  0.271(U2  4sec  6),JUJFA  (3-19) 

The  F  -insolation  at  level  4  is  therefore 
A 

FA(4)=  FA[1  -  .271(U0  .sec  0)*303]  (3-20) 

A      A  2,4 
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Of  the  energy  F.(4),  the  fraction  1  -  R  is  transmitted  through  the 

cloud  top  which  interfaces  at  level  4.   Hence  the  downward  transmitted 

F  -insolation  through  the  cloud  top  is 
A 

FA(4)i  =  FA(1  -  Rc)[l-.271(U2j4sec  6)*JUJ]         (3-21) 

Within  the  cloud  layer,  the  water-vapor  path  in  layer  (4,6)  is  subject 
to  the  mean  diffuse-path  augmentation  factor  5/3  (after  Katayama,  1966). 
Thus  the  energy-absorption  A(4,6)  is  given  by 

A(4,6)  =  FA(4)i  *.271[|(U4j6  +~)]*JUJ  (3-22) 


In  (3-22),  the  water-vapor  parameter  \]     is  the  water»vapor  mass  equivalent 
of  the  cloud-droplet  absorption.   The  determination  of  the  parameter  Ur 
which  is  spread  vertically  over  the  cloud-layer  (4,8)  is  modeled  after 
estimates  of  cloud-droplet  absorptivity  a_T   [given  by  Korb  et  al.(1956)] 
and  makes  use  of  the  modelling  assumption 

aCLD  "  -271  {[!<U4,8  +  V1'303-  'I  U4,8''303}  <3"23> 

The  method  of  computation  of  U_/U.  Q  is  detailed  in  Appendix  C  as  a  function 
of  U,  8  for  the  latitude  bands  0-40,  40-65,  65-90. 

With  U-  known  at  each  gridpoint,  the  following  three  computations  of 
insolation  are  made 

A(4,6)  =  FA(4)i  {.271[|(U4>6  +/)]*303}  (3-24) 

A(6,10)=  FA(4)i  {.271[|(U4)1Q+  Uc)  ]  *3°3-[f  (U4>6+  ^l'303}       (3-25) 

14   -  FA(4)i  {l-0.271[f(U4  lfl  +  Uc)]'303}  (3-26) 
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These  three  computations  represent  respectively  (i)  the  absorbed  insola- 
tion in  the  layer  (4,6),  (ii)  the  absorption  in  layer  (6,10),  and  (iii) 
the  transmission  to  the  earth's  surface.   The  last  of  these  is  subject 
to  the  multiplicative  earth-absorption  factor  P  (of  form  similar  to 
Eq.  3-8) 

P  =  (1  -  ag)/(l  -  ag  Rc)  (3-27) 

This  factor  P  when  multiplied  by  I   J  of  (3-26)  gives  the  F  -penetration 
of  the  earth's  surface  under  cloudy  skies. 

In  summary,  the  gridpoint  absorption  of  F  -energy  in  earth  now  becomes 

A 

the  cloud-weighted  sum 


IA(10)  =  (CL)(P  of  3-27)  (I  c  of  3-26)+(l  -  CL)(IAQof  3-18)    (3-28) 
Similar  cloud-weighted  sums  of  absorption  of  F  -energy  should  be  formed 

A 

for  the  layers  (6,10)  and  (2,6).   In  the  latter  connection,  it  should  be 
recalled  that  the  F  cloudy-sky  result  in  layer  (2,6)  is  made  up  of  the 

A 

two  contributions  (3-19),  (3-22),  the  former  of  which  has  no  cloud-reflec- 
tivity factor  (1  -  R  )  since  it  pertains  to  A(2,4)  only. 
3(d)   Summary  of  solar-heating  effects 

The  results  of  Sections  3b,  3c  are  now  combined  to  get  a  resultant 
insolational-absorption  at  earth.  This  is  done  by  summing  the  two  contri- 
butions Ic(10)  of  Eq.  3-12  and  1.(10)  of  (3-28). 

o  A 

Within  the  atmospheric  layers  (2,6),  (6,10),  only  the  F  wavelengths 

contribute  to  absorption  and  therefore  to  temperature-change  through 

2 
°T^o  cs   _  A(2,6)  cals/cm  /hr  /"»_?q\ 

St(2'6)  "     iooo ^-^) 
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J&6.10)    =  A(6,10)    cal/cm2/hr 

g   p 

The  first  of  these  temperature-changes  involves  the  cloud-weighted  com- 
posite absorption  rate  formed  from  A(2,6)  of  (3-15)  and  from  A(2,4), 
A(4-6)  of  Eqs  (3-19,22).   Similar  weighted  means  are  needed  for  A(6,10). 
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4.   Formulation  of  net  cooling  rates  by  terrestrial  radiation 

(a)  Clear-sky  areas 

It  should  be  recalled  that  Table  1  contains  the  format  of  a  "radiative" 
sounding  having  input  data  defined  at  the  levels  listed  in  Table  3. 

Table  3.   Levels  and  data  used  in  the  "radiative" 
soundings. 


Levels 

10 

8 

6 

4 

2 

1 

0 

Temperature 

ho'  Ts 

T8 

T6 

T4 

T2 

Tl 

To 

W.V.mass 

U10=° 

U8 

U6 

U4 

U2 

Ul 

uo 

C0_  mass 

C10=° 

C8 

C6 

C4 

C2 

Cl 

co 

The  value  of  water-vapor  reduced  absorber-mass  U_,  (k  =  5,....,0)  was  defined 

in  Section  2(d),  as  was  also  the  corresponding  value  of  C   .   These  values 

have  been  incremented  upward  relative  to  the  ground  as  indicated  by  Eqs. 

2-12,  2-11,  respectively.   In  these  equations  the  notation  U(2k)  has  been 

treated  interchangeably  with  the  more  compact  notation  U  ,  used  here,  and 

in  the  remainder  of  this  section. 

Consider  the  radiative  temperature  change  in  air.   By  the  first  law  of 

thermodynamics,  this  is  given  by 

* 

(4-1) 


•k  ~k 

5T     (F6  "  F10> 


St 


o.4n  1222  c 


g    P 

for  the  layer  (6,10).   For  the  layer  (2,6),  the  corresponding  formula  is 

given  by 

*    * 
(F„  -  F,) 

(4-2) 


dT 

*   * 

(F2  "  V 

5t  _ 

0.4H  10°°  C 

g        P 
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The  numerators  of  (4-1),  (4-2)  which  are  contained  within  parentheses 
are  understood  to  be  flux-divergences  from  the  atmospheric  columns 
(6,10)  and  (2,6)  over  the  gridpoint  having  the  indicated  "radiative" 
sounding. 

Virtually  all  of  radiation  packages  in  the  presently  used  large- 
scale  circulation  models  in  the  United  States  [cf.,  Manabe  and  Strickler 
(1964),  Oliger  et  al.  (1970),  Arakawa  et  al.  (1969)]  make  use  of  Yama- 
moto's  (1952)  water  vapor  flux-emissivities  as  a  basis  for  their  radia- 
tive flux  calculations.  These  absorptivities  e(U,T)  as  defined  below 
are  virtually  independent  of  temperature  in  the  normal  range  of  tempera- 
tures (300K  >  T  ^  220K),  but  at  temperatures  T  <  220°K  the  wave-averaged 
flux-emissivity  function  defined  in  (4-3)  after  Yamamoto  (1952)  and 

others,  has  the  form    °° 

[  dBv 

J  tl-T  (L  U)]  — i  (T=220K)dV  (4-3) 

6(U,T  <  220)  =  ^ H^2 dT 


dB 
dT 


)     -^  (V,T)  dV 


It  is  seen  that  the  numerator  of  (4-3)  has  its  temperature-dependent  terms 

3 
fixed  at  220K,  while  the  denominator  equals  4St*T  even  when  T  <  220K. 

Thus  (4-3)  introduces  a  temperature-dependent  flux-emissivity  for  tempera- 
tures T  <  220K,  reflecting  the  effect  of  increased  absorptivity  by  the 
rotational  band  of  atmospheric  water-vapor  which  takes  effect  in  view  of 
Wien's  law  at  these  low  temperatures.   In  the  discussion  which  follows, 
the  absorptivity  will  be  restricted  to  be  due  to  water-vapor  alone,  until 
the  subject  of  C0„  is  specifically  introduced.   Although  analytic  functions 
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have  been  determined  for  6(U,T)  for  both  water  vapor  and  C02,  it  is  useful 
to  examine  schematic  versions  of  the  Yamamoto  method  of  net  flux  calcula- 
tions through  key  levels  of  the  sounding  (Table  3).   Schematic  net  flux 
depictions  are  shown  in  Fig.  3  for  the  levels  10.,  6.,  2.  in  cases  (A), 
(B),  (C),  respectively. 
4a.  (A)  Net  flux  at  level  10. 

Sasamori  (1968)  has  empirically  recomputed  6(U)  after  Yamamoto 's  (1952) 
values  for  the  temperature- independent  part  of  atmospheric  radiational 
soundings,  schematic  versions  of  which  are  shown  in  Fig.  3.  The  tempera- 
ture-independent emissivity  6(U),  after  Sasamori  (1968)  is 

e(U)  =  .240  log1Q(U  +  .010)  +  .622  (4-4) 

in  the  temperature- independent  region  T  ^  220K.   For  T  <  220K,  Sasamori 
gives  for  the  temperature-dependent  emissivity 

€(U,T)  =  (8.34T-3531og10U-44)  D-034551og10U-.705  (4.5) 

Here  U  represents  the  final  U-value  relative  to  the  reference  level 
considered  to  be  level  10,  so  that 

i(U,T)  =  (8.34T-3531og10Uo-44)  Uo--03455log10Uo-75         (4-6) 

Corresponding  to  this  value  of  U  ,  a  temperature-dependent  e-isopleth 
converges  into  the  lower  left  corner  of  each  panel  of  Fig.  3,  i.e.,  6  -»  1.0 
as  B(T=0)=  0,  where  B  =  Sc*T  is  the  Stefan-Boltzmann  flux. 

The  trapezoidal  summation-rule  gives  a  representation  of  the  hatched 
(downward  flux)  area  of  Fig.  3A  as 

Fj(10)  =  .5  |e(8-10)(B10-Bg)  +  [6(8-10)+  6(6-10)] (Bg-  B6) 

+  [e(6-l0)+  e(4-lO)](B6-B4)+[e(4-iO)+  e(2-i0)] (b4-b2) 

B(Tj) 
+  [6(1-10)+  6(2-10)]  (B^H  2^  6(0-10, T)dBJ  (4-7) 

B=6 
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where  £  in  the  final  term  is  the  function  of  U,T  listed  in  (4-5).   The 

4 
notation  6(8-10)  denotes  e(U0-Uir.)  etc,  while  Bir.  =  St*T,  _  , 

-9      -2   -1  o  -4 
St  =  5.67332  x  10  cal  cm  "  hr    K   and  so  on  for  all  of  the  other  black- 
body  fluxes  (BR,  B-.,  B,,  B  ,  B. ,  B_),  distributed  along  the  sounding  of 

Table  3. 

r  - 

Note  that  the  final  integral  \  6  dB  of  (4-6)  may  be  evaluated  by  the 

integration 

T=T 

\  €  dB  =  8.34  UQ  io-.034551og1()U0-.705  ^   T.3531ogloV.44  4St*T3dT 

0 

.3531ogU  -.44 
-  W1{8.34U0-^531ogU0-.705.   >  I 


10  0 


It  finally  follows  that 


,Bi 

\     '*  dB  -  St*Tl4  *  1{V0>V   *  .3531og^oV3.56 
o 


(4-8) 


with  e(UQ,T1)  defined  in  (4-5). 
Hence  F.Q  is  given  by 

ho  '  Bio  -  F.(10)  (4"9) 

with  the  latter  term  given  by  Eq.  4-6,  4-7. 
B.   Net  flux  at  level  6 

In  this  case,  the  desired  net  flux  is  represented  by  B(T1Q)  minus  the 
slant-hatched  area  above  the  sounding  curves  of  Fig.  3B.   Thus  the  result 
may  be  modeled  after  the  results  4-6,  and  4-7 
F*  =  B   -.5  {e(6-8)(B8-B6)+[e(6-8)+  e(6-10)] (B1Q-  Bg) 

+  e(4-6)(B6-B4)+[e(4-6)+  e(2-6)](B4-  B2) 

+[e(2-6)+  6(1-6)] (B2-   Bx) 

+   (St>TA   'e<U0,6>Tl  >   *  .353   JgU^^  3.56|  (4"10) 
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Here  €(U,T)  is  defined  in  (4-5)  while  e(  |u2k  -  U&  |)  is  defined  in  (4-4). 

C.  Net  flux  at  level  2 

Fig.  3C  is  relevant  for  this  computation.  The  net  flux  at  level  2 
may  therefore  be  written  as  follows  by  analog  with  Eqs.  4-8,  4-6,  4-5, 
4-4: 

F2*=  St*T1Q4  -  .5  |e(2-4)(B4-  B2)  +  [€(2-4)  +  6(2-6)] (Bg-  B4> 

+  [£(2-6)  +  6(2-8)]  (Bg-  B6)  +  [€(2-8)  +  €(2-10)  ]  (B^-  Bg) 

+  €(102)(B2-B1)  +  St*^4*  -€(U0>2Tl)*  <353  lo^>2+  3<56J     (4-1D 

4 
where  B,  =  St*T 

■k  *  * 

With  F_  ,  F,  ,  Firj   formulated  in  terms  of  radiative-sounding  parameters 
by  Eqs.  4-7  and  4-9,  4-10  and  4-11  respectively,  layer-mean  temperature-change 

rates  in  K(hr)  may  now  be  calculated.   Before  this  is  done,  however,  it  is 

*    *     * 
convenient  to  introduce  the  modifications  in  F«  ,  F,  ,  Fin  resulting  from  the 

inclusion  of  CO  flux- emission. 

* 

D.  Inclusion  of  C0_  absorber-mass  in  F  calculations 

CUN 
The  "area"  integrals  \     €(U,T)  dB  depicted  in  Fig.  3  have  been  computed 

U=0 
by  integrating  along  a  two-path  sounding  curve  away  from  the  reference  level 

(Fig.  3A,B,C).   When  absorption  is  due  both  to  water  vapor  and  C0_,  the 

product-law  of  transmissivities  must  be  employed.   In  doing  this  the  trans - 

missivity  of  the  two  absorbers  acting  jointly  is  expressed  as 

TF(U,C)  =  Tp(U)  TF(C) 

Then  the   flux-integral  becomes 

,UN'CN  >V 


Yi  =  )     [i-TF(u)TF(c)]dB  =  ^  jU-yu)]  +  Tp(U)   eF(c)j  dB  (14-12) 

U,C=0  U=0,C=0     "w.v.flux"  "CO, 


II/-TI     M 
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The  first  term  on  the  right  is  the  water-vapor  flux  computed  using  for 
example  Eq.  4-7  and  ignoring  the  emission  by  CCL.   The  second-term, 
denoted  "CO  ",  represents  the  effective  emission  due  to  CO  not  already 
counted  in  the  water-vapor  flux  computation.   The  added  emission  is 
expressed  schematically  on  the  Yamamoto  chart  by  an  increase  of  emissivity 
at  all  points  of  the  radiative  soundings  (except  at  the  reference  level 
points).   Thus  at  any  temperature  T,  the  required  joint-emissivity  is 
obtained  by  replacing  that  due  to  water-vapor  already  computed,  €  (U) , 


by 


Swc5Vu)  +Aewc(u'c>- 


The  last  term  includes  the  15  |0m  band  emitted  flux  of  C0„  which  has  not 
already  been  counted  in  emission  by  the  rotational  band  of  water  vapor. 
Sasamori  (1968)  has  computed  this  15  \Jm   band  "overlap"  region  emissivity, 
and  developed  the  empirical  correction  shown  in  (4-13,  4-14)  below: 

Aew,C(U2k>C2k>  =  C  «W  *  VC2k>  (4"13) 

Aew(,(2k)  =  .07262  {l-. 62556  (U2k+. 286)  *26j  |log  C2k+  1.064J      (4-14) 

Thus  in  considering  the  flux-modification  for  the  inclusion  of  CO,,, 
one  simply  computes  the  added  emissivity-values  Ae   of  (4-14)  at  each 
sounding  point.   In  the  summation  Eq.  4-7,  for  example,  where  €(8-10) 
appeared,  this  term  is  simply  replaced  by  the  function 

VU8  "  V  +  A£WC[<V  <V>  <V  V1 
where  the  function  A€   term  uses  the  simultaneous  values  Ug  -  U1Q  and 

CQ  -  C._  stored  from  the  radiative  sounding.  The  15  pm  effective  emissivity 
8    10 

Ae   of  Eq  4-14  is  considered  to  be  temperature-independent  (Sasamori,  1968) 
WC 
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in  the  range  T  ^  T-^220K,  in  which  range  the  water-vapor  emissivity  € 
was  temperature-averaged.   Hence  in  Eq,  4-7,  the  final  term  of  F*(10) 
becomes 

St*Tl4'VVTl)(.3531og^V3.S6>  +  6WC  of  (4"W)1  (4"15) 

where  A€  _  in  the  last  term  of  the  bracket  in  (4-15)  involves  substitu- 
WC 

tion  of  UQ,  CQ  for  U,C  in  (4-14). 

*    *     * 
With  the  revised  formulations  for  Fo  »  fa  '  F10  ifc  is  aPProPriate 

to  compute  the  flux  divergences  in  the  layers  (2,6),  (6,10)  and  finally 

Fio*  ■  Bio  -  Fwc<10>4  (4"16) 

which  gives  the  terrestrial  radiative  loss  at  the  surface  with  clear  skies 
above. 

In  addition,  only  a  simple  programming  alteration  is  required  to 

*     * 
include  the  C0_  modifications  to  recover  the  flux  divergences  F,  -  F10 

*    * 
and  F_  -  F,  necessary  for  the  computation  of  the  layer  temperature- 
2.         o 

change  rates,  Eqs.  4-1,  4-2. 

* 
4(b)  Modification  of  net-flux  F  with  cloudy  skies. 

As  before,  the  cloud-cover  amount  CL  is  given  by  Eq.  2-9.   Section  4(a) 

has  dealt  with  the  clear-sky  fraction  (1-CL)  over  gridpoint  (i,j),  and 

the  remaining  fraction  CL  therefore  represents  an  overcast  (extending 

vertically  from  a  =  .8  to  a  =  .4)  over  the  remainder  of  the  grid  area. 

*    *    * 
As  before,  we  are  interested  in  F«n  ,  Ffi  ,  F«   (but  with  overcast)  over 

the  gridpoint.   The  relevant  schematic  flux  diagrams,  after  Yamamoto 

(1952),  are  depicted  in  Fig.  4.   In  Fig.  4,  all  emissivities  are  based  upon 

the  function  e   as  computed  using  Eqs.  4-4  and  4-14  in  proper  combination 

of  absorber  masses  relative  to  the  reference  level. 
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Fig.  4.   Schematic  area-depictions  of  net  flux  at  reference  levels  10 
and  2,  when  an  overcast  cloud  extends  between  levels  8  and  4  in  the 


vertical.  .The  horizontally-hatched  area  is  F-._  ; 
area  is  F   .   The  net  flux  F,   is  zero,  since  it 
body  cloud.   F  *  is  numerically  equal  to  F.n  . 

o  1U 


the  slant-hatched 
lies  within  the  black- 
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4(b) A.   Net  flux  F-„   at  earth,  with  overcast  clouds. 

The  result  for  this  case  is  based  upon  Fig.  4    and  is  given  by  the 
horizontally-hatched  area: 

Fio*  -  <Bio-  V[1  -  -5  Wio'W  (4"17) 

This  is  a  cooling-rate  at  the  earth-surface. 

At  the  same  time  the  net  flux  at  level  8  (=  the  overcast-base)  is  an 
identically  equal  magnitude,  so  that  the  layer  8-10  is  in  radiative  equilibrium. 
However,  the  cloud  base  receives  net  flux  warming  exactly  equal  to 


* 


F1Q  ,  i.e., 


■k  -k 

Fg  =  -  F1Q   ,  a  warming  (4-18) 


4(b)B.   Net  flux  at  level  6. 

Since  a  layer  of  approximately  50  meters  of  average  cloud-drop  content 
may  be  regarded  as  a  black  body  (Brunt,  1939)  at  the  mean  temperature  of 
the  cloud  layer,  it  follows  that  level  6  is  characterized  by 

F6*  =  0  (4-19) 

4(b)C.   Net  flux  at  level  2. 

The  design  of  the  result  in  this  case  is  suggested  by  Fig.  4.    Thus 

* 
F_  becomes 

F2*   -  V    .5|V<2-4><VB2)+  V^XVV   + 
4.0e(u        T  ) 
2V.353  lo8;oP02+3.56  +  Aewc%,2  C0,2>  °f   <4"14>]}  <4"20> 


after  application  of  Eqs.  4-11,  4-12  to  the  truncated  area  of  Fig.  4. 
This  net  flux  is  a  much  smaller  value  than  that  of  Fig.  3C  because  of  the 
cold  black  body  interposed  at  level  4/. The  reduced  net  flux  due  to  clouds 
is  an  effect  often  referred  to  as  "cloud-contamination"  in  problems  of 
remote  sensing  of  the  atmosphere. 
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4(c)   Terrestrial  radiative  cooling  rates  in  key  atmospheric  layers 
(i)   Upper  layer  (2,6) 


The  radiative  flux  divergence  in  the  cloud-covered  case  is 
r*   of  4-20)  -  F6 


(F  *  of  4-20)  -  F  *(=  0)  (4-21) 


This  result,  appropriately  cloud -weighted  with  the  clear  sky  result 

*     * 
F0   -  F,   of  Section  4(a),  then  gives  the  terrestrial  radiative  cooling 

*     * 
in  the  layer  (2,6).   The  cloud-weighted  F«  -  Ffi  may  then  be  used  in 

connection  with  Eq.  4-2  to  yield  a  temperature-change  rate  due  to 

terrestrial  radiation.   The  latter  may  then  be  combined  with  the  insola- 

tional  absorption  of  3-29  to  obtain  an  overall  radiative  change  in  the 

layer. 

(ii)   Lower  layer  (6,10) 

The  corresponding  result  for  this  cloud-covered  layer  is 

(F(*.  0)  -  Fg*  -  Fg*  -  -  F1Q*  (4-22) 

where  a  negative  result  indicates  a  net  flux-convergence,  or  a  warming 
rate.   The  result  of  (4-22)  covers  the  layer  (6,8).   To  this  may  be  added 
the  zero-change  in  layer  (8,10),  so  that 

F6*  "  F10*  "  "  F10*   of  (4"17)  (4'23) 

(iii)   The  interface  level  a  »  1.0 

* 
The  terrestrial  radiative  transfer  across  the  earth-surface  is  F.q 

which  must  now  be  cloud-weighted  in  the  usual  manner,  i.e. 

F1Q*(i,j)  =  CL[F1Q*  of  (4-17)]  +  (1-CL)  [T  lQ*   of  (4-16)]  (4-24) 
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5.   Maintenance  of  the  heat  budget 
5(a)   In  the  atmosphere 

At  each  of  the  odd  levels  in  Fig.  1,  the  diabatic  heat-sources  of  the 
radiative  type  discussed  in  Sections  3  and  4  are  included  in  the  thermo- 
dynamic equation 

|f(ncpTk)  +  v  •  (ttcpVV  +  S3  (TV>A>  "  v'\ 

Here  Q,  is  the  diabatic  heating  (cooling)  rate  per  gram  at  level  k  includ- 
ing that  due  to  radiation.   Q,  should  also  include  the  latent  heat  released 
at  level  k  at  time-step  t.   The  subject  of  latent- heat  release  and  more 
generally  of  the  atmospheric  water  vapor  balance  is  not  included  in  this 
paper. 
5(b)  At  the  earth's  surface 

The  earth's  surface  is  considered  to  be  a  heat  reservoir  of  zero  heat 
capacity  over  land,  and  of  infinite  heat  capacity  (a  perfect  insulator) 
over  the  ocean.  The  case  of  a  solid  ice  sheet  may  be  considered  to  be  a 
predictable  compromise  between  the  land  versus  ocean  cases.  These  earth- 
surface  vertical  heat-flow  requirements  lead  to  the  following  surface 
boundary  conditions,  essentially  after  Langlois  and  Kwok  (1969)  henceforth 
abbreviated  to  LK. 
5(b) A.   Over  land 

The  external  radiative  surplus  (deficit)  at  level  10  is 

* 
I   -  F 
10    10 

With  the  land  surface  considered  to  be  of  zero  heat  capacity,  this  surplus 

must  be  diffused  upward  in  the  form  of  sensible-heat  conduction  T   and 

latent  heat  of  evaporation  E  at  the  earth  surface,  across  land: 

Iio  '  Fio*  =  r  +  E  =  r[i+i]  (5-i) 
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In  (5-1),  r  =  r/E  is  the  so-called  Bowen  ratio,  given  empirically  by  Sellers 
as 

r  =  9.6(sin  Cp)2  -  7.93  (sin  cp)  +  2.0  (5-2) 

sin  Cp^  0.4  for  all  Cp 

In  this  treatment,  a  direct  formulation  of  E  over  land  is  not  achieved; 
rather  only  an  indirect  computation  making  use  of  E  =  T/r  is  employed. 

The  determination  of  the  sensible-heat  flux  T   is  based  upon  its  large- 
scale  parameterization  [(after  Deardorff  1966),  LK(1969)],  and  a  simultaneous 
micrometeorological  parameterization.   Thus  LK  give  the  large-scale  para- 
meterization as 

r  =  -PioCpKf    ^ 

r  -  PioV      a*(99-  es)>k  -  ^  (5"3) 


1  + 


Z9   "   Z10 


and   the  micrometeorological   equivalent   as 

r  "  P10CpVs(Tg   -  V  (5"4) 

In  (5-4),  V_  is  related  to  the  surface  geostrophic  wind  through  V  =  - 8 |V   | 
S  o      go 

+2.2  in  mps.  In  (5-3)  K  ,  a  and  y  are  turbulent-transfer  coefficients 
whose  values  have  recently  been  documented  by  Kaitala  (1972).  The  value 
of  0   is  related  to  the  estimated  T   =  Tg  by 

es=Tio<lSoo>-286  <5"5) 

where  T  -  T  is  the  actual  temperature  lapse  (compatible  with  Eq. 5-3)  within 
g   x 

the  constant  flux  boundary  layer. 

A(i)  Overland  boundary  adjustment  of  surface  temperature 

The  combination  of  Eqs  5-3,  5-4  leads  to  the  solution  for  Tx  at  the 
top  of  the  boundary  layer: 
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( y  y  tcDvsig  -  k,c]  +  k  e9 

In  (5-6)  6  is  related  to  Tq  by  a  formula  of  form  (5-5).   If  one  starts  over 
land  with  T  =  T„,  a  solution  for  T  emerges  and  a  value  of  T  is  available 
from  (5-4).   Application  of  the  Bowen-ratio  technique  then  provides  a  value 
of  E.   With  first  iterates  of  T   and  E  available  for  use  in  (5-1),  an  adjust- 
ment in  surface  temperature  T. _  is  achieved  as  follows 

Fio*<y  •  ho  -  ^ +  b  (5-7) 

F10*<V  ■  B10  -  Fd  •  F!0*<V  +  <i>Ts(Tg  -  TS>  <5-8> 

If  one  makes  use  of  Eqs  5-7,  5-8  for  correction  of  T  in  5-6,  a  revised 
r  of  (5-5)  results  and  sequentially  there  is  a  second  iterate  T  ^  from  5-7, 

5-8.   This  process  must  converge  rapidly  in  view  of  the  large  coefficient 

3        3 
4  T   or  4  T  ,  ,  depending  upon  the  step  in  the  iteration  procedure.   When 

b       8^- 
the  several  iterative  steps  are  complete,  a  heat  balance  exists  at  the 

earth's  (land)  surface,  which  has  realistically  determined  the  surface  tem- 

*       * 
perature.   Induced  changes  in  F,   and  F_  have  resulted  from  the  recalcula- 

b       Z 

tion  of  the  earth-surface  temperature  to  be  T  (instead  of  T  ).  These  become, 

8  5 

respectively,    in  the  case  of  clear-skies 

/tf6*  =    [B(T      )-  B(TS)]    [1  -    6^(6-10)]  (5-9) 

AF  *  =    [B(T      )    -   B(T_)]    [1   -   €      (2-10)]  (5-10) 

z  gn  b  WC 

where  n  is  the  iterate-step  and  e   is  given  by  Eqs,  (4-4,  4-14) 

With  overcast  skies  only  the  minor  correction  below  applicable  in  the  layer 

(8,10)  occurs 

AF*(8,10)  =  [B(T  )  -  B(Tg)]  [1  -  ewc(UMQ)]         (5-11) 
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A(ii)   Over  the  ice-free  ocean. 

Here  it  is  assumed  that  Eqs  5-5,  5-6  still  apply,  but  that  T  is  known 
from  the  sea-surface  analysis,  which  is  held  constant  during  the  prediction 
process.   In  this  case,  a  unique  value  of  T  and  therefore  of  T  results. 
It  is  no  longer  necessary  to  make  the  Bowen  ratio  assumption  for  the  surface- 
layer  evaporative  heat-transfer  rate.   Instead  the  analog  of  (5-3)  is  formed 


E  =  fio  .8.-  9.  T~^Z  (5"12> 


K* 

qx  "   q9 

e  -  e 

*      9         S 

1  +  a  (z  -z    > 

*9      10 

Z9   "Z10 

and  is  coupled  with  the  analog  of  (5-4) 

E  =  ^oVs^sg"  %?  (5"13) 

Eqs.  5-12,  5-13  yield  the  result 

a  .  *»,  +  <V  Z!0>  CDVsg  (5.14) 

"x    K  +  (Z9-  Z10)  CD  Vs 

with 

,._   (AE  -  BE/TK) 
q  =   -622  e 2_  (5-15) 

Sg  n  -  e(AE"BE/Tg) 

AE  =  21.656;   BE  =  5418.   The  last  relationship  is  the  Clausius-Clapeyron 
equation. 

In  this  case  (5-1)  may  be  rewritten 

Fio*  "  ho  -  <E  +  r'  (5"16) 

If  F    of  (5-16)  differs  from  that    computed  using  Eq. (4-16),  the  differ- 
ence I  -  Flfl  should  be  used  to  increment  both  E  and  T  in  the  relative 
proportion  of  the  magnitudes  of  the  latter  pair. 
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A(iii)  Over  ice-covered  ocean 

In  this  case,  the  surface-to-air  boundary  condition  may  be  stated  as 

r  ■  ^o  ■  Fio*  -  b<Tw  -  y  <5-i7> 

-2   -1   -1 
where  b  =  ice-conduction  coefficient  (697.8  erg  cm  sec  deg   ) 

T  =  271. 2K  is  the  freezing  point  of  sea  water 

w  -or 

T  =    the  desired  temperature  at  the  ice-air  interface 

The  evaporative-heat  transfer  E  has  been  omitted  from  5-17.   However,  Eqs. 

5-3,  5-4  are  still  valid  and  it  is  still  necessary  to  obtain  T  in  order  to 

x 

evaluate  T  by  5-4.   Accepting  T  =  T  as  the  first  iterate,  it  is  possible 
to  adjust  F-n  (T  )  of  (5-17)  as  was  done  in  Eq.  5-8.   If  only  a  single 
iteration  is  required,  an  analytic  formula  for  T  results  [Kaitala  (1972); 
LKQ969)],  namely: 

£°  +  FT7-  '  &C 
9   10 
T  =  — - (5-18) 

X    1 (1000  .286  Q_0 

Z  -Z.  C.9TT  }  +   P10  p  DVSU  Q) 

9     b 


where 


D  "  ho  +  bTw  +  4St*TS4  "  F10*(TS) 

P  -  piocp  .ce9.  es) 

[1  +  a  - — = — J 
I  =   Cp/CCp  +  b  +  4St*Tg3) 


°r '  PioWs  <5-19> 
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With  T  known  in  terms  of  T  ,  one  may  iterate  successively  upon  the 

values  of  the  latter  using  5-8  and  obtain  convergent  values  of  T  as 

8 

well  as  of  the  sensible  heat  transport  T.      With  T  replaced  by  a  more 

S 

realistic  value,  T  ,  it  is  relevant  again  to  replace 

o 

F6*by     F6*(Tg) 
F2*by     F2*(Tg) 

Fio*by  Fio*(y 

following  the  method  of  Eqs.    5-9,    5-10. 
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APPENDIX  A 
Vertical  extrapolation  of  the  mixing  ratio  profile. 

When  London's  (1957)  climatological  soundings  included  mixing  ratios 
q(p)  to  levels  of  p  =  0.3TT,  the  interpolated  values  of  q (a)  (where  10.0a  = 
9.,  7.,  5.,  3.)  were  subjected  to  a  best-fit  against  pressure  in  the  loga- 
rithmic form 

ln  aM  =  x    ln  p_  .  x  ln  5  (A„1} 

q5       p5     .5 

Eq.  (A-l)  was  tested  for  eight  latitude  belts  on  the  Northern  Hemisphere  and 
over  four  seasonal  mean- soundings  at  each  belt.  The  parameter  X  is  related 
to  the  correlation  coefficient  r(Y,X)  (where  Y  =  ln  q/q,-  and  X  *  ln  p/p,-), 
which  was  well  in  excess  of  r  ^  0.99  in  each  sounding.   This  indicated  an 
effective  method  of  extrapolation  in  the  form 

2- =  (^-)X  (A-2) 

q5     P5 
as  also  suggested  by  Smith  (1966).   This  relationship  permitted  extrapolation 
of  individual  moisture  soundings  to  p  ■  .3TT  when  q(p)  values  existed  at 
p  =  .9TT,  .7tt,  .5TT  at  analysis  points. 

When  moisture  values  q(p)  exist  also  at  p  ■  .lTT,  a  relationship  of  form 
(A-l)  was  also  found  to  give  close  verification  at  O  ■  O.lTT  in  terms  of  a 
best-fit  power  law 

-2-  =  (^-)X  (A-3) 

q3      p3 
with  X  determined  using  the  four  pairs  of  (q/q~>  p/Po)  values  from  the 
analysis  covering  each  gridpoint. 
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APPENDIX  B 
Solar  undepleted  insolation  in  spherical  coordinates. 

The  incoming  undepleted  insolation  is  given  by 

F  =  S[sin  Cp  sin  6  +  cos  Cp  cos  6  cos  h]  (—)2  (B-l) 

m 

where  (r/r  )  is  listed  by  Julian  date  in  the  Smithsonian  Tables  (List,  1963) 
m 

number  169.   Likewise,  the  declination  6  is  listed  in  that  same  source,  but 
may  also  be  computed  from 

sin  6  =  .39785  sin[4. 88578  +  .0172D  +  .03342  sin(.0172D) 
-.001388  cos(.0172D)  +  .000348  sin(.0344D)  -.000028  cos  (.0344D)]  (B-2) 
The  hour  angle  h  on  day  D  is  given  in  terms  of  GMT  (denoted  by  Z)  by 
means  of  the  relationships 

t  =  M0D[(Z  +  2Tl724)  +  24'  24] 

cos  h(Z,X)  =  cos  [|£  (t-6)  -  H]  (B-3) 

The  expression  in  brackets  in  (B-l)  is  cos  6,  where  6  is  solar  zenith 
angle  at  latitude  Cp,  6  and  hour  angle  h,  the  latter  given  by  B-3.   If 
cos  6  of  (B-l)  becomes  negative,  it  is  set  to  zero. 
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APPENDIX  C 

Determination  of  the  absorber- equivalent  water  vapor 
mass  U  for  droplet-absorption  in  clouds 

In  clouds,  the  water  droplet-absorption  has  been  modeled  using  the  dif- 
fuse-path absorptivity  of  Manabe-Moller  (cf.,  a__   of  3-23).   In  this  paper 
the  value  of  a_T_,  due  solely  to  cloud  droplets  has  been  taken  after 
theoretical  computations  of  Korb,  Michalowski  and  Moller  (1956),  and  sub- 
sequently used  also  by  Manabe  and  Strickler  (1964).  These  values  are  as 
listed  below  in  Table  4. 


Table  4.   Absorptivities  (a   )  of  cloud  droplets 

CLD 


oud-type 

aCLD 

R 
c 

Equivalent 

Low 

.035 

.21 

0-39N 

Middle 

.020 

.48 

40-64N 

High 

.005 

.69 

65-90N 

Since  only  a  single  cloud-layer  has  been  invoked  in  the  radiative  model 
presented  here,  the  larger  absorptivities  associated  with  low  clouds  (lar- 
gest droplets)  have  been  associated  with  an  equivalent  latitude  band  0-39N; 
while  intermediate  values  of  aPT   have  been  identified  with  an  intermediate 
latitude  band,  40-64N.  Finally,  the  minimum  value  a„Tn  has  been  ascribed 
to  the  latitude  band  65-90N.   The  values  ar   of  Table  4  are  due  solely  to 
absorption  by  either  the  liquid-water  or  the  solid-ice  content,  depending 
upon  the  cloud  type  considered.  The  values  aPT   of  Table  4  are  then  to  be 
considered  absorption-means  for  the  latitude  belts  under  consideration. 
No  specific  use  has  been  made  of  the  values  of  R  ,  although  it  seems  reason- 
able to  ascribe  these  reflectivity  values  to  the  same  latitude  bands  for 
which  a_T  -values  are  proposed  in  Table  4. 
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The  absorptlvities  of  Table  4  apply  over  the  full  energy  range  of  the 
solar  spectrum.  If  cloud-water  absorptivities  are  to  be  modeled  as  water- 
vapor  equivalents  occurring  solely  in  the  F. -wavelengths,  then  each  value 
of  aPTn  must  be  multiplied  by  (1./.349)  on  the  left  side  of  (3-23} .  Thus 
one  obtains  the  water  vapor  equivalent  of  the  cloud-water  absorptivity  in 
the  form 

a*  -  2.8653aCLD  =  {[fo^g  +  Uc)]-303-  [fu4>8r3°3}(.271)     (0-1) 

using  a  Manabe-Moller  parameterization  similar  to  (3-23).   In  (C-l),  the 
absorption  in  clouds  has  been  spread  over  the  full  cloud-depth  which  con- 
tains U.  0,  a  known  absorber  mass  of  water  vapor  as  determined  from  the 
4,8 

sounding.   Eq.  C-l  may  be  solved  for  U  /U,  g  using  the  appropriate  value 

■k 

of  a     =  2.8653aPT      and  Table  4.      The  desired  solution  assumes   the   form 

(1  +!c__)-303  =  1+JL  (_3Z5_r303 

U/    «  L        .271    \  R   ; 

4,8  4,8 

which  becomes 

4,8  4,8 

where  the  exponent  1^0.303  =  3.3003 

A  graphical  solution  of  (C-2)  for  U£  as  a  function  of  U4  g  is  presented 
in  Table  5  as  a  varies  over  the  low,  middle  and  high  cloud  (i.e.  latitude) 
types  of  Table  4. 
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%* 

uc(i) 

0.01 

0.1 41 70 

0.0?. 

0.18936 

0,05 

0  . 2  6  c  8  6 

0.10 

0.4C852 

0.15 

0.50529 

0.20 

0.53995 

0.30 

o.  i i 7 at? 

0.50 

0.9  8617 

0.  80 

1.25740 

i.CO 

1.48119 

1.20 

1.05219 

1.40 

1.  3  11-32 

1.60 

1.96648 

1.80 

2.11306 

2.00 

2.25403 

2.  20 

2.39016 

2.40 

2.52206 

2.60 

2.65022 

2.80 

2.77499 

3.00 

2 .89675 

Uc(2) 

j  . 

05118 

0. 

0  7  252 

0. 

12003 

0. 

17^51 

0. 

22673 

0. 

26  939 

0  . 

344  73 

0. 

4  7267 

63  54  8 

0. 

7322^ 

0. 

82271 

0. 

908  22 

0. 

9(<^74 

1. 

,06793 

1. 

14331 

1. 

,2162  3 

1 

.2  3609 

1, 

, j 55 34 

1 

,t2297 

1. 

,  4  3  o  5  4 

U   (3) 
c 


0.00740 
0.31155 
0.02103 
0  .0^329 
0.0436^ 
0.05294 
0.06955 
0. 058  2  6 
0.135  20 
0.15735 
0.17822 
0. 15  8  JO 
0.21691 
0.235  10 
0.  2t>2o7 
0.269  70 
0.28626 
0.30240 
0.31816 
0. 33357 


Table  5.     Water-vapor  equivalent  U     for  droplet  absorption  in 
relation  to  the  pressure-scaled  water-vapor  absorber  mass  U,    fi 
(in  the  sounding  layer   (4,8)).     Case  U   (1)   corresponds  to        ' 
low-cloud  water-drop  absorptivity;  '  U   (2)   to  mid-level 

cloud-drop  absorptivity;   U   (3)   to  high-clSud  absorptivity. 
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